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We compare the prediction of Skyrme-like effective Lagrangians with data for electromagnetic 
form factors of proton and consider the possibility of fixing the parameters of these higher-order La- 
grangians. Our results indicate that one or two-parameter models can lead to better agreement with 
the data but more accurate determination of the effective Lagragian faces theoretical uncertainties. 
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I. INTRODUCTION 

The Skyrme model Q , despite is relative successes Q , can only be considered as a prototype of an effective theory 
of QCD. Indeed, large N c analysis 0, Q suggests that bosonization of QCD would most likely involve an infinite 
number of mesons. If this is the case, then taking the appropriate decoupling limits (or large mass limit) for higher 
spin mesons leads to an all-orders Lagrangian in derivatives of pion fields. For now however, fixing the form of the 
effective Lagrangian from an exact low-energy limit of QCD seems to be out of our reach and the alternative has 
been to propose simple effective Lagrangians ja> IE an d rei Y 011 a f ew nucleons properties to set the parameters they 
depend on. From that point of view, a rather stringent test for such models lies in whether they could accurately 
describe the data for the electromagnetic form factors of nucleons for moderate values of momentum transfer. 

The electromagnetic form factors of the semiclassically quantized SU(2) skyrmion were studied systematically by 
Braaten et al |8| but this first attempt did not take into account relativistic kinematical corrections which are important 
for momentum transfers Q 2 > 1 GeV 2 . These correction were implemented to the original Skyrme model predictions 
by Holzwarth H G3 using the prescription of Ji Actually, Holzwarth also introduced a second correction to 

the Skyrme model to account for vector meson effects either by introducing a multiplicative factor to reflect the 
contributions of the poles coming from these mesons, or by addin g th e vector mesons as dynamical degrees of freedom 
in the Lagrangian. Yet, in spite of remarks and suggestions in jloL Il2l Il3j. no analysis regarding how higher-order 
Lagrangians could reproduce the nucleon electromagnetic form factors or conversely, how the form factors could help 
construct a more accurate effective Lagragian. 

In this work, we calculate the electromagnetic form factors for a class of higher-order (Skyrme-like) models intro- 
duced in . For simplicity, we limit our numerical analysis to one and two-parameter models with a proper treatment 
of relavistic effects, and evaluate to what extent the experimental data of the electromagnetic form factors of proton 
can suggest a form of higher-order Lagrangians or discriminate among possible candidates. Note that neutron data 
are ignored here since they are plagued with large uncertainties. 



II. THE SKYRME MODEL 



Let us first introduce the Lagrangian density for the Skyrme model 

C s = -^Tt(L^) + ^Tt[L^L u ] 2 (1) 

where is the pion decay constant, L M is the left-handed chiral courant = U^d^U and the Skyrme constant e is a 
dimensionless constant. U is a SU (2) field related to the pion field tt by U = exp(2«7r • t/F^). The field configurations 
with finite energy must satisfy the boundary condition 

U(r, i) — ► 1 for|r|^oo. (2) 
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These configurations fall into topological sectors characterized by 



1 



B= ^ / d J xdet{L?} 



£l jk 
~48tt 2 



(3) 



a topological invariant taking integral values. 

Skyrme interpreted this topological invariant as the baryon number. Accordingly, the lowest-energy B = 1 sector 
is identified with the nucleon. In this sector, the lowest-energy field configuration is given by the hedgehog ansatz 



U(r) = exp [it ■ rF(r 



(4) 



where F(r) satisfies the boundary conditions F(0) = it and F(oo) = 0. 

With the convenient change of scale, we can use 2y/2/eF 7r and F n /2-^2 as units of length and energy respectively 
and rewrite the Lagrangian density Q as 
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1 



1/1 



(5) 



where f^ v = [L^LJ\. A pion mass term 14] is usually added to account for the chiral symmetry breaking observed 
in nature 



r 



m 2 F 2 

(Till - 2) 



(6) 



This term serves as a regulator for the magnetic radii and form factors (see (|36H ) of nucleon otherwise they would 
diverge [T^ . 

Using the hedgehog ansatz , one obtain the mass of the static skyrmion 



M = - J d 3 x£ s 



drr 2 < F 



sin 2 F sin 2 F 
2 — + — 
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sin 2 F 



2/3 2 (l-cosF) 



(7) 



where r has now been rescaled and /3 = 2v / 2w 7r / 'bF„ . The stable static soliton is obtained by minimizing the mass 
and requires to solve the chiral equation 



(l + o) 



F 1 sc 
F" + 2— - 2 — 



,c sc F' 
- + --2 — 

s r z r 



-/5 2 s = 0. 



(8) 



with the boundary conditions F(0) = ir and F(oo) = for 5 = 1. For simplicity here, we used a = sin 2 F/r 2 , 
s = sinF, c = cosF. 

Fluctuations around this static soliton soliton should be quantized. Quantization of the skyrmion is usually per- 
formed with the introduction of spin and isospin rotation matrix as a collective coordinate 0. The spin/isospin 
rotations of skyrmion takes the form 



U{r,t) =Ai(t)U(r)A(t) 
with A(t) an arbitrary time-dependent SU(2) matrix. Substituting JJjJ in JSJ), one gets 



L = -M + ITi[d t A^d t A} = -M + 7 ^ ^ 



where M is defined in (J7| and 



J: 



8tt / 2\/2 
T I ~^F, 



2 \ f 

- \ I r 2 dr sinF(r) 



sinF(r) F , 2 



(9) 



(10) 



(11) 



is the moment of inertia of the skyrmion. Here I is the spin or isospin of the nucleon. The parameters F v and e are 
fixed using two experimental input (mass of the nucleons or else) . One is then able to reproduce the static properties 
of baryons within a 30 % accuracy [3] ■ These methods can also be generalized to extensions of the Skyrme model as 
we will show in the following section. 
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III. ALL-ORDERS SKYRMIONS 



The Skyrme model is more a prototype for an low-energy pion interactions than a full effective held theory. Higher- 
order terms are expected to appear in addition to JSJ and Q but in its most general form, the Lagrangian would involve 
a increasing number of terms at each order in pion held derivatives making the treatment practically intractable. One 
of us has proposed a special class of models whose energy density, assuming the hedgehog ansatz, is at most linear in 
F' 2 |(|. This requirement is sufficient to determine a unique term to each order in derivatives and turns out to have 
deeper geometrical meaning 

The static energy density coming from the Lagrangian of order 2m in derivatives of the held takes the form 



a" 1 " 1 [3a + m(6 - a)] 



(12) 



where a = sin 2 F/r 2 and b = F' 2 . Using the hedgehog ansatz, the hrst two terms arise from the non-linear a and the 
Skyrme terms 
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(13) 



while the third term leads to 



£ 2 = -£ 2 = -IsrtfafV = a [a + 2b] 
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"•£3 — T^Tr/^,,/^ f x 



(14) 
(15) 



as for the term proposed by Jackson et al 0] to allow for the dynamics of the u) meson in the Skyrme model. 

Generalizing to all order, the static energy associated to this class of all-order Lagrangian can be written in a simple 
form 



£ = 2J h m £ m = 3x(a) + (b- a)x'(a) 



(16) 



where a specific model is characterized by a choice of the parameters h m or equivalently of the function x( a ) = 

00 

^2 h m a m and x'( a ) — iF- Yet, x( x ) is no t completely arbitrary. Requiring that a unique soliton solution exists sets 

m— 1 _ 

some constraints on x{ x ) : 

d_ 

dx' 
d ( x{xY 



-X(x) > 
< 



dx V x 



d ( 1 d , . , 



for x > 0. 

The mass of the soliton, including the pion mass term, is 



M = Air 



( F « 



\2V2~t 



OO 

^ J r 2 dr{3 X (a) + (b-a) X ' {a) + 2/3 2 {l - cos F)} 



(17) 



and leads to the generalized chiral equation 



X (a) 



F' sc 
F" + 2 2 — 



+ ax" (a) 



,nC SC F' 

F' 2 - + — - 2— 



P 2 s = 0. 



(18) 



In that context, the Skyrme Lagrangian corresponds to x( a ) = a + \- The moment of inertia of the soliton also take 
a simple form 



_ 8tt ( 2V2 




3 V e 3 F, 



r 4 dra[2x'(a) + (b-a)x"(a)}. 



(19) 
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Most all-orders model depend on more than two parameters. In the next section, we analyze the behavior of the 
electromagnetic form factors for a few models in the hope that these could help fix the h m coefficients and obtain a 
better agreement with the the experimental data in general. 

IV. ELECTROMAGNETIC FORM FACTORS 

The electromagnetic form factors of the proton in the Breit frame, for spacelike momentum transfer q 2 > 0, are the 
Fourier transforms of its electric charge and magnetic moment densities: 

oo 

G e ( -I 2 ) = \ J dr {Bo i^+B^r)} jo (qr), (20) 



CW V) = M„/*{l^+ IBlW }^ (21) 



with 

B (r) = — sin 2 (F)F', (22) 

Bx (r) = g (^j u*a [2 X > (a) + (b a) X " (a)] (23) 

are the baryon density and moment of inertia density respectively. Here j n is the spherical Bessel function of order n 
and Mjv is the nucleons mass. 

Both densities (|22() and (|23[1 are normalized 

oo oo 

J drB (r) = J drB x {r) =1 (24) 
o o 

while the electromagnetic form factors satisfy the normalization condition 

(0) = 1, G P M (0) = v P (25) 

where 



and 



3 \2I 

is the magnetic moment of the proton with baryonic square radius 



e 2 F 2 

o 



drr 2 B {r). (27) 



However, the definitions l|2(J|) and (|21|l only hold in the Breit frame moving at velocity v with respect to the nucleoli 
rest frame where the chiral profile F(u) is computed. A correction for this Lorentz boost must be applied. Ji 01 has 
proposed a simple prescription to circumvent this difficulty: 



G E tf) = G%(£). (L'S) 



G M (q 2 )= 1 - 2 GZ (4)- <- >! » 
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where G]f and G^J are given respectively by (|20() and i|21|) and 7 is the Lorentz factor 

^(l- 2 )- 1 -!^ (30) 

with the nucleon mass M . 

Unfortunately, the boost transformations l|28|) and (|29[) violate the so-called superconvergence rule 

q 2 G E ,Ai{q 2 ) -> 0, for </ 2 -> 00. (31) 

which is expected to hold for electromagnetic form factors. Indeed the limit q 2 — > 00 in the Breit frame corresponds 
to q 2 = AM 2 in the rest frame and generally G^ M (4M 2 ) does not vanish. 

A possible approach to restore superconvergence is to relax the condition that M must take the value of the 
nucleon mass and instead allow M to vary in (|30() in order to get the best agreement with the data of G p m /{p p Gd) 
at the highest available values of q 2 01 • However for the models under study here, this procedure turned out to be 
unsatisfactory. The value of M ensuring superconvergence caused the ratio of the electromagnetic form factors to 
become too suppressed in the large-g 2 limit. On the other hand, the models which we propose as candidates for the 
description of QCD at low energy are not expected to hold for large values of q 2 . Since the superconvergence rule 
seems to be too restrictive, it will not be applied here. We will instead promote M as parameter and adjust its value 
to provide a better fit of high-g 2 data. 



V. RESULTS AND DISCUSSION 



The parameters of the Skyrme model, e and F n , are usually set with either one of the following two methods: 
(i) adjust e and F v to obtain the mass of the nucleons (939 MeV) and of the A resonance (1232 MeV) or (ii) set 
according to its experimental value (186 MeV) and adjust e to reproduce the nucleon-A mass split (295 MeV). 
For comparison purposes, we adopt the second method also used in [jj, IToj and assume that the pion mass takes its 
physical value. 

We consider here two simple extensions of the Skyrme model that fall into the class of models described in Section 

eh 

a 2 

Model A: Xx{a) = a + — + ca 3 (32) 

3 

Model B: Wo) = a + , Cl ° , (33) 

(1 + cia) 

Model A adds a single term (and a single parameter c) of order six in derivative of the field to the Skyrme model 
whereas Model B which depends on two additionnal parameter has a rational form which hopefully could allow to 
reproduce poles due to vector mesons. This latter model is a generalization of a model introduced by Jackson et al 
Q with ci = i and 02 = f ■ 

Fixing the model parameters e, c, c\ and ci as well as the scale parameter M is a tedious procedure which requires 
a few steps: First, we choose a set of parameters (c for Model A or c\ and C2 for Model B) and solve the differential 
equation ijTSjl for massive pions and use the method described above (inputs are F^ and nucleon-A mass split) to fit 
for the appropriate value of e. In the second step, we compute the form factors, compare them with data and adjust 
M to minimize the x 2 of the ratio of the electromagnetic form factors. The whole procedure is repeated with different 
sets of parameters until we get the configuration e, c, M or e, Ci, C2, M with lowest \ 2 ■ 

The results for the electromagnetic form factor are presented in Fig. ^We get e = 5.03, c ~ ^ and M — 1.20 GeV 
for Model A and e = 3.40, c\ ~ ^, C2 — 1 and M = 1.20 GeV for Model B. The results for the Skyrme model are also 
shown for comparison (e = 4.25 and M = 1.66 GeV). According to Fig. ^ both models present improvements over 
the Skyrme model, especially for the ratio Ge/h p G v m of the form factors. Despite its simplicity, Model A seems to 
overcome Model B. However both models exhibits a rapid divergence of the magnetic form factor of proton. 

In view of these results, we must conclude that although the models examined here provide clear improvements, 
they cannot reproduce the data for the magnetic form factor adequately. The problem originate to their inability 
to mimic the pole coming from the vector mesons. Computing these pre dictions for an arbitrary number of Skyrme 
model extensions being prohibitive, we resort to another approach |9Lll0j|. The p— meson effects are incorporated by 
multiplying the form factors l|20|l and l|21|l by 
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FIG. 1: Electromagnetic form factors of proton from Model A: c 
(dashed lines) and Skyrme model: e = 4.25 (dotted lines), 
squares) . 



e = 5.03(full lines), Model B: ci 



The data are from Q,[l9, 20] (full diamonds) and [l3 E^il (full 



1, e = 3.40 



with M p = 770 MeV. The electromagnetic mean square radii of proton then become 



6 dG P E (-q 2 ) 



E G(0) dq 2 
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q i=o M 2 e 2 F 2 



drr 2 {B {r)+B 1 (r)}, 



(35) 
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2 n 



dr r 



4r 2 B (r) 
e 2 F 2 l 



■lB x {r) 



(36) 



This correction introduced an additional parameter whose allowed values go from the purely pionic model (A = 0) to 
the mesonic dominance (A = 1). It should be noted that the parameter A does not have any effect on the ratio of the 
electromagnetic form factors. A similar effect could also be included to take into account the 10— meson effects, but 
corresponding factor would involve a second A parameter and our calculations indicates that this additionnal factor 
does not leads to noticeable amelioration of the behavior. 

Figure |21 presents the electromagnetic form factors as obtained from Model A and B with pole effects introduced 
through (|34|l . All the models parameters as well as the parameters M and now A are fitted using the procedure 
described above to minimize the \ 2 of the magnetic form factor. The best accord with data are acheived for values of 
the parameters e = 4.74 c ~ ^ with M — 1.31 GeV and A = 0.73 for Model A and e = 3.40, c\ ~ | and ci ~ 1 with 
M = 1.26 GeV and A = 0.72 for Model B. The results for the Skyrme model are found to be consistent with those 
of refs. 0,0] with M — 1.42 GeV and A = 0.75. Clearly, the inclusion of vector mesons effects produce a significant 
improvement G v M j [l v Gb over the results of Fig. ^ Again both extensions of the Skyrme model, Model A and B, are 
better at reproducing the data with Model A being slightly superior despite its relative simplicity. Furthermore, as 
shown in table [I] the addition of higher orders terms have little effect on other static properties of the proton with 
these choice of parameters (Note that the magnetic moment of proton predicted by the models are all lower than 
its experimental value but one has to remember that these properties are subject to quantum corrections 23]). The 
values of A indicate a strong dominance of the p— meson as in [9|, [lfj . 
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FIG. 2: Electromagnetic form factors of proton from Model A: c = e = 4.74 (full lines), Model B: ci = 
(dashed lines) and Skyrme model: e = 4.25 (dotted lines). The data are the same as for Fig..0 



c 2 = 1, e = 3.40 





Skyrme 


Model A 


Model B 


Exp. 


M (GeV) 


1.43 


1.31 


1.26 




A 


0.75 


0.73 


0.72 




Hp (mat) 


1.78 


1.79 


1.80 


2.79 


r% (fm) 


0.684 


0.701 


0.706 


0.870 ± 0.008 


r P M (fm) 


0.826 


0.841 


0.843 


0.858 ± 0.056 



TABLE I: Static nucleon properties and parameters for models of FIG|5| The experimental data come from [24112 



Our results therefore indicates that the data are best reproduced by a higher-order effective Lagrangian in the low- 
energy limit of QCD . Unfortunately, as we have noticed, the determination of a more acurate effective Lagrangian 
faces theoretical uncertainties which remain to be addressed. First, the boost prescription (|28|l and l|29|) . not being 
compatible with superconvergence, prevents us from reproducing results for q 2 > 10(GeV/c) 2 . An adequate boost 
prescription holding account of superconvergence should allow to extend the analysis to higher momentum tranfer. 
Secondly, we have only explored the possibilities of two types of Skyrme model extensions. These could not mimick 
the vector meson effects satisfactorily without the artifact of equation (|34|l . Finding a form of Skyrme- like lagrangian 
that would allow to bypass the approach in (|34[) and avoid the introduction of the A parameter remains a challenge. 
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